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A nonlinear two-fluid stochastic hydrodynamical descrip- 
tion of velocity and concentration fluctuations in sedimenting 
suspensions is constructed, and analyzed using self-consistent 
(SC) and renormalization group (RG) methods. The advec- 
tion of particles by velocity fluctuations is shown to be rel- 
evant in all dimensions d < 6 . Both RG and SC analyses 
predict a strong reduction in the dependence of velocity fluc- 
tuations on system-size L relative to the L 1//2 obtained in 
the linearized theory of Caflisch and Luke [Phys. Fluids 28, 
785 (1985)]. This is an important step towards resolving a 
ten-year old puzzle in the field. 

Sedimentation Q is at once a rich and complex phe- 
nomenon in colloid science and a frontier problem in 
nonequilibrium statistical mechanics. The average sedi- 
mentation speed w aod is determined by balancing the driv- 
ing force (gravity) against the dissipative force (viscous 
drag). Despite the long-ranged (1/r) nature of the hydro- 
dynamic interaction between the colloids, an adequate 
mean-field theory (J] of v scd is available for low concen- 
trations. There is, however, no satisfactory theory of the 
fluctuations of the velocity and concentration fields in a 
steadily settling many-particle suspension, even at zero 
Reynolds number. Experiments show clear nonther- 
mal fluctuations in the particle motion, with the char- 
acter of random stirring on a scale substantially larger 
than the particle size. The fluctuations Sv in the settling 
speed are of the same order as the mean v BCd . Caflisch and 
Luke (CL) || showed in a linearized hydrodynamic anal- 
ysis that, for a sedimenting system of linear dimension 
L, the assumption of purely random local concentration 
fluctuations led to long-ranged velocity fluctuations with 
Sv ~ L 1 / 2 . Experiments, however, find no dependence 
of Sv on L In this letter we assume a source of 

randomness at small scales, the origin of which we com- 
ment on briefly below, and examine its consequences at 
large scales. In detail, we use a dynamical renormaliza- 
tion group (DRG), general scaling arguments, as well as 
perturbative self-consistent methods to analyze the ad- 
vection of concentration fluctuations in steady-state sed- 
imentation by the velocity fluctuations which they them- 
selves produce. The DRG method (see Ma and Mazenko 
jd|), was applied to the closely related problem of the 
randomly stirred fluid by Forster, Nelson, and Stephen 
(FNS) Q . Scaling and self-consistent approaches give a 
family of possible results, the simplest and most phys- 
ically appealing of which accounts for the experimen- 
tally observed lack of size-dependence. The DRG results, 



which are to first order in e = 6 — d where d is the dimen- 
sion of space, present a clear trend towards suppressing 
the size dependence, but simply setting e = 3 in these 
first-order calculations leaves a residual size dependence, 
albeit much weaker than the L 1 / 2 of CL H. 

We coarse-grain the dynamics on scales comparable to 
the random-stirring scale l s . This obliges us to use a 
stochastic description with a noise source and an effec- 
tive bare diffusivity or hydrodynamic dispersion coeffi- 
cient p3[ Do ~ v scd £ s . We then calculate quantities such 
as the velocity fluctuation on scales 3> i s using the follow- 
ing two-fluid hydrodynamic description of sedimentation 
along z. Our coupled equations for fluctuations v and 
c in the velocity and concentration fields respectively, at 
zero Reynolds number, are 

Br 

— + Av Vc = D V 2 c + 6(r,t) (1) 
at 

r/V 2 Vi = m R g8cP iz (2) 

Here the pressure field has been eliminated by impos- 
ing incompressibility via the transverse projector Py = 
Sij — ViVj (V 2 ) -1 . Equation [J is the advection-diffusion 
equation for a scalar field with A = 1 . The bare collective 
diffusivity D and the Gaussian noise term 9, with zero 
mean and 

(6(r,t)0(O,0)) = N V 2 S(r)S(t) (3) 

have the same origin — in general a combination of ther- 
mal fluctuations and hydrodynamic dispersion. If the 
noise is of purely chaotic origin, as in Ladd's athermal 
simulation then N and D are parameters which 
must be determined either experimentally or by the anal- 
ysis of numerical simulations. For purely thermal noise, 
the ratio N /D is fixed by the fluctuation dissipation 
theorem (FD) Q while D a ~ k B T/i]a with n the sol- 
vent viscosity and a the radius of the colloids. Eq. ||, 
which expresses the balance between the driving by grav- 
ity and dissipation by viscosity, describes how concen- 
tration fluctuations produce velocity fluctuations with 
mug the buoyancy-reduced colloidal weight. In Eqs. |l| 
other possible nonlinear terms, e.g. those arising from 
concentration-dependent mobilities, are readily shown to 
be subdominant relative to the advective nonlinearity 
V.vc @. 
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If we neglect the nonlinear term (i.e. set A = ), 
we can solve the equations to find the result of |J: 

Q 



then 



(4) 



The nonlinear term causes these induced velocity fluctua- 
tions to advect further the concentration field. Since flow 
advection is well known to distort and suppress fluctua- 
tions (l^] , the ultimate size-dependence should be weaker 
than that predicted by Eq. |[ 
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FIG. 1. The renormalization group flow below the up- 
per critical dimension in the coupling constant (A), noise 
anisotropy (Fjv) plane. The IR stable fixed point which vio- 
lates the fluctuation - dissipation theorem (F^A*) is repre- 
sented by the black dot at the intersection of the dashed lines 
representing the fixed point values of A and Tjv. 

Let us summarize first the simple scaling and self- 
consistent arguments. Note first that within the above 
linear treatment the dynamic exponent relating spatial 
and temporal rescaling is z = 2 and the structure factor 



for concentration fluctuations S(k) =< 



Ck 



>^ for 



wavenumber k — > 0. The nonlinearity V.cv, whose ad- 
vective form guarantees that there will be no graphical 
corrections to it under renormalization, alters these con- 
clusions dramatically. Rescaling (isotropically for sim- 
plicity) r — > br, t — » b z t, c — > b Xc c, v — > b Xv v, normaliz- 
ing to retain the form of the first term in ([[]) , and insisting 
that the coefficient of the nonlinearity remain unchanged, 
yields Xc + z + ^- = Xv+z — 1 = 0. We can identify the 
upper critical dimension d c above which the nonlinearity 
is irrelevant by noting that for d > d c , linearized be- 
haviour should obtain: Xc = —d/2, z — 2. This gives 
d c = 6, which is some indication of why a reliable answer 
in three dimensions is going to be hard to find. The same 
results can obtained from a one-loop self-consistent cal- 
culation of the renormalized diffusivity D(k,uj) and noise 
spectrum N(k,u>) jlE|. In either case, if we assume that 
there is no singular renormalization (except possibly log- 
arithms) of the noise spectrum 111], we get, for d < 6, 



z = d/3, and S(k) ~ 1/D{k) ~ k 2 ^ 3 . In three dimen- 
sions this is equivalent to the result S(k) ~ k of [ jl9| . 
While this last is but one of a family of possible results 
that arise out of the scaling analysis, it has a particu- 
lar physical appeal, which we discuss at the end of this 
Letter. 

We turn now to the DRG calculation (lljjl^] of the 
renormalization of the dynamical parameters N, D, and 
A. The RG recursion relations to one loop order near the 
upper critical dimension d = 6 are 



1(T 4 A 2 (£) {1.0 + 9.4rAr(£) - 3T%(£)} 
= 2.62 x 10~ 6 A 2 (£) x 



1 

dr N 

di 

[1 - 5.6T N (£) - 163. 8F^) - 23ir^( 

— = ~10- 4 N(£)X 2 (£) x 
di 

[0.55 + 5.5T N (£) - 5.5T 2 N (£)] 

^ = -6.72 x 10- 5 X 2 (£)T D (£) x 
[1 + 9.63Fjv(£)] 



(5) 

(6) 
(7) 
(8) 



We have allowed explicitly for different noise strengths 
N z and N± for wavevectors along and normal to z re- 
spectively. Here, the dimensionless coupling constant is 

A = * /2 , and e = d — 6, while Tn = N±/N z and 

Tf) = D±/D z are the anisotropy factors of respectively 
the noise and diffusion constants. We selected the dy- 
namical exponent such that D z has trivial scaling. Figs, 
and show the resulting RG flow lines for d > 6 and d < 6 
respectively. For d > 6, the nonlinear coupling constant 
X(£) scales to zero for large £. Linearized hydrodynam- 
ics is applicable and (v 2 )l oc L A ~ d + const in agreement 
with Eq. ||. For d < 6, on the other hand, we find a sta- 
ble fixed point (A*,r^). The noise is highly anisotropic 
at the fixed point (r^ w 0.0608), while A* oc sfe. The 
noise amplitude N(£) scales to zero according to Eq. [?], 
as does the diffusion constant anisotropy ( Eq. ||). This 
means that the ratio N{£)/D{£) of the noise to the diffu- 
sion constant scales to zero at the fixed point signalling 
the absence of an FD theorem at large lengthscales be- 
low d = 6, irrespective of the origin of the noise. There is 
a FDT respecting fixed point which is unstable to any 
bare deviation of the noise from its FDT determined 
value. The linear theory, which in effect assumes an 
FD theorem, thus overestimates the magnitude of con- 
centration fluctuations at large lengthscales. This flow- 
induced noise-suppression is highly anisotropic. This is 
already clear from the fact that is a small number, 
so concentration fluctuations with wavevectors along the 
sedimentation direction are much larger than concentra- 
tion fluctuations perpendicular to the sedimentation di- 
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rection. A similar conclusion can be drawn from the dy- 
namical scaling exponents. The relaxation rate ui{q z ) for 
concentration fluctuations along the z direction depends 
on the wavevector q z as ui(q z ) oc q z , with z\\ the cor- 
responding scaling exponent, while in the perpendicular 
direction w(q_i_) oc |q±| Zi . If we use e = 3 in Eq. [| we 
find zu = 0.91 while z± = 1.93. In-plane relaxation of 
concentration fluctuations is thus nearly purely diffusive 
while out-of-plane relaxation is much faster than diffu- 
sive. 
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FIG. 2. The renormalization group flow above the up- 
per critical dimension in the coupling constant (A), noise 
anisotropy (Fjv) plane. The IR stable fixed point is that de- 
scribing the free theory with a finite shift in the anisotropy 
parameter Tjv = N±/N z . This fixed point model obeys the 
fluctuation - dissipation theorem. 

To make quantitative predictions we use the dynamical 
scaling relations 



D z ,x(q z ,±,u) = q z *±G z ,x(u/q z z Y) 



where a 



-2 + zii + 7, f3 z 



z\\, Pi 



(9) 
(10) 
-2 + z± 



with 7 = 0.83 for e = 3. The scaling functions F Zt ±(x) 
and G Zt ±(x) are all finite at x = but for large x, they 
behave as power laws: 



limF z ,±(x) cx ie q/z "^ 

a;— > 

limG z j_(x) cx x^.J-Aii.-l 
The full velocity fluctuation spectrum is then 
<k 2 (q^)l> = 

N z (q,Lu)q 2 z +N ± (q,Lu)qi 
lo 2 + [D z {a,w)ql + D^io)qlf 
|vi(q,«)|) = 

q\q\ \ N z {q,uj)ql + N^u)ql 

oj 2 + (D z {d,u>)q 2 z +D x (c l ,oj)q 2 ± f 



(11) 
(12) 



A- 



A 



(13) 



(14) 



with A = (mug/rj) . Integrating over q and u> while 
using the dynamical scaling relations given by Eqs. |s|JTo| 
gives our final result for the size-dependence of the ve- 
locity fluctuations in d = 3. 



\v z \ 2 ) L L 1+ ^- z «V2-t 



Vj_ ) L OC 



L l- 7 



(15) 
(16) 



The trend towards a substantial reduction of the size- 
dependence is gratifying. If we simply extrapolate our 
first-order in e results to d = 3 we find (\v z \ 2 )l oc L 6 
and (|v^| 2 )i oc L 017 . Comparing Eqs. [l5] and |l6| with 
experiment we first note that velocity fluctuations di- 
verge considerably more weakly than linearly along the 
z direction and hardly at all along the perpendicular di- 
rection so Eqs. |l5] and O are closer to the observed size 
dependence than Eqs. ||and ||. While experiments re- 
port no size dependence at all, our theory really is only 
reliable near d = 6 and a naive extrapolation of a one 
loop DRG calculation to d — 3 is questionable. An- 
other interesting method of comparison between exper- 
iment and the present theory is in the anisotropy fac- 
tor (|vjj 2 )/(|uj| 2 ). In the linearized theory this quantity 
is equal to 1/3 (Eq. ^) while the non- linear theory pre- 
dicts that this ratio vanishes as L~ QA3 with system size 
(Eqs. [l5|, |l^). Ladd's simulation gives 0.12 for this ratio 
and recent experiments 0.27 and 0.26 |||l7]]. All of these 
measurements thus suggest that the ratio is indeed less 
than 1/3. 

Before closing, let us note that the naive self-consistent 
results in d = 3 can also be obtained by dimensional argu- 
ments reminiscent of those used in the theory of homoge- 
neous isotropic turbulence [ pp[ , but with the assumption 
of a scale-independent velocity rather than dissipation. 
z = 1 or equivalently D(k) ~ 1/k means that the dif- 
fusivity at a length scale k~ 1 is given by that length 
times a constant which must on dimensional grounds 
have the units of velocity. In the absence of external, 
thermal noise, the only candidate for this constant is 
w sod itself. Making the reasonable guess that all physi- 
cal quantities, in the Pe — -> 00, Re — ► limit, must then 
be built from ?; acd and wavenumber leads immediately to 
D(k) ~ w BCd /fc, Sv ~ v sod and the mean-square displace- 
ment of particles ~ v 2 od t 2 upto logarithmic factors and 
volume-fraction dependent coefficients. 

Finally, we summarize the achievements of this work. 
First, starting from the correct, coarse-grained, nonlin- 
ear stochastic hydrodynamic equations, which in princi- 
ple are a complete description of the problem, we have 
shown clearly that the origin of the disagreement between 
theory || and experiments [^|J^] lies in the neglect of an 
important nonlinearity, viz., the advection of concentra- 
tion fluctuations by the fluctuations in the flow. Our 
work thus emphasizes that both the renormalized sedi- 
mentation speed |Q and fluctuations therein should be 
calculated in parallel with, and not subject to assump- 
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tions about, the structure factor, and shows that stochas- 
tic models are a way of doing this. Second, our iden- 
tification of d = 6 as the upper critical dimension for 
the problem shows that the effects of this nonlinearity 
are likely to be (a) very large and (b) rather difficult 
to calculate analytically in three dimensions. Third, our 
dynamic renormalization group calculations in an expan- 
sion in e = 6 — d indicate a clear trend in the direction of 
reduced system-size dependence. Finally, with one addi- 
tional assumption, we can use self-consistent and dimen- 
sional arguments which give results which imply at most 
a logarithmic size-dependence, and make further predic- 
tions for experiment. A reliable calculation directly in 
three dimensions remains the main challenge for theory. 
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